Torque Theory of Anisotropic Superconductors with No 
Phenomenological Parameter in determining Vortex Core Size 

Daichi Kubota^, Nobuhiko Hayashi^'*-^, and Takekazu Ishida^'*"'^ 
^Department of Physics and Electronics, 
Osaka Prefecture University, 1-1, Gakuen-cho, 
Naka-ku, Sakai, Osaka 599-8531, Japan 
^ Nanoscience and Nanotechnology Research Center (N2RC), 
Osaka Prefecture University, 1-2 Gakuen-cho, 
Naka-ku, Sakai, Osaka 599-8570, Japan 
^ CREST(JST), Sanban-cho, Chiyoda-ku, Tokyo 102-0075, Japan and 
^Institute for Nanofabrication Research, 
Osaka Prefecture University, 1-1, Gakuen-cho, 
Naka-ku, Sakai, Osaka 599-8531, Japan 
(Dated: January 13, 2013) 

Abstract 

The contribution of vortex core has been taken into account properly in constructing a torque 
theory for multiband superconductors. We employ the prescription of describing internal magnetic 
field in the vortex lattice by Hao et al. and by Yaouanc et al. to derive a torque formula as a natural 
extension of a preceding London theory. In marked contrast with the preceding model, our novel 
formula does not contain a phenomenological parameter r], which prevents us from obtaining a true 
upper critical field Hc2 by analyzing an experimental torque curve. The parameter r] was originally 
introduced to take care of the uncertainty in determining the vortex core size Furthermore, 
we reveal that the rj value is universally scaled by anisotropy 7, magnetic field and Hc2 due to 
field dependence of ^y. This may revitalize the single-band Kogan model in combination with a 
universal function r]{-y, B, Hc2) instead of a constant r/. 
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Thanks to extensive studies on the superconductivity of a magnesium diboride and iron 
arsenides, the torque as the angular derivative of a free energy has been recognized as a 
powerful tool for investigating superconducting anisotropy. For example, MgB2 can be char- 
acterized by isotropic three-dimensional tt band and anisotropic two-dimensional a band. 
Discovery of high-Tc cuprates strongly suggested that lower dimensionality is a key prereq- 
uisite for inducing high-Tc superconductivity. A recent study by Yuan et al. [l| is indicative 
of an opposite view because they found a rather isotropic superconducting anisotropy in 
(Ba,K)Fe2As2. As a sensing probe of anisotropy, the torque would play an important role 
for revealing the nature of various new superconductors. 

Theoretically, the magnetic torque of an anisotropic superconductor has been investigated 
for a long time. Originating from the local-limit model of the torque, it evolved into a model, 
which takes care of the nonlocal effect for extending to a variety of superconductors. In the 
early London model of the torque ^ , a vortex core is not considered seriously but is simply 
assumed to confine at the region shorter than a coherence length ^. The later Kogan model 
[sl discusses not only the anisotropy of a penetration depth A but also that of ^ independently, 
but it still remains to be a local model. Yaouanc et al. ^ computed the Fourier components 
of the magnetic field in a high-K type-II superconductor containing an ideal vortex lattice 
for analyzing the local magnetic field in /iSR and neutron diffraction measurements. Their 
method revealed the effect of the vortex core for determining local magnetic field. Brandt 
developed an analytical method using full Ginzburg-Landau calculations for ideal vortex 
lattice with any spatial symmetry [a]. The method developed by Yaouanc et al. gave a 
useful too. iu aua.yz,ug the vortex state iu NbSe, single erysta. J Q and the effeet of the 
multiband nature on the local field profile in /iSR experiments [8| . 

There is an issue to be amended in the London model because the cutoff in the reciprocal 
Fourier space (G-space) at 27r^~^ coming from the vortex core is introduced rather a priori. 
There is still room for improving the London theory since it is recognized that the core 
size (or its G-space equivalent) depends on the magnetic field. In this Letter, we attempt to 
derive the torque formula on the basis of the discrete Fourier- analysis method of the internal 



magnetic field developed by Yaouanc et al. 



Before explaining our new theory, we first summarize the development of the preceding 
torque theories. The torque is directly obtained by differentiating a free energy with respect 
to angle 9 between crystalline axis and the magnetic field. The free energy density F is 
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given by [9| 



87r$o 



(1) 



where hJG), as in Ref. 



9|, is the component of the magnetic flux in the G space, B is 
the magnetic flux density, and $o is the quantized flux. We take coordinates with z— axis 
parallel to magnetic field. The magnetic flux hz{G) is derived by 

<l>o(l + A..G2) 
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where Ajj is a penetration depth tensor, Ajj = X^niij (Aq = A^m^, Ac = A^mc) is represented 
by rrixx = XHij. cos^ 9+mc sin^ 9, m^z = ('ria — mc) sin 6' cos 6', rriyy = rria, and rrizz = sin^ + 



rrir cos^ ^ 



io|. 



Suppose that the relation Hd ^ i/ ^ i/c2 is satisfied in the intermediate field so that 
an effective vortex spacing L can be expressed as ^ ^ L ^ A. Using the dimensionless 
reciprocal lattice vectors g = LG with L = ($o/-B)^^^, we approximate free energy density 
within the zeroth and first-order expansion with respect to g. Actually the free energy 
density is approximated as |lO| . 

B^m^.L"^^' 1 
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where sums up g (g 7^ 0). A free energy density is given from Eq. (|3]) as 10| 
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where a phenomenological parameter 77 (~ 1) was first introduced 10| and has been anxious 
about the clarification of its physical meaning. Thus, one finds a single-band torque formula 
of anisotropic superconductors as {2] 
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where e{9) = (sin^ 9 + -f"^ cos^ 9)^/"^. 

Multi-band effect in anisotropic superconductors leads to the separation of the degener- 
ated electron, an.sotropy parameters rnto two. i.e.. an.sotropy i,r eohererrce length = 
and anisotropy in penetration depth 7^ = \c/^a [21, l3|. The magnetic field in G space is 
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assumed to be zero at G > Gmax so as to fulfill the cut-off of the vortex core at ~ ^ in real 
space. One expresses the free energy density of Eq. (|3]) as 

^ %Bmzz /■^■""'^ dG^dGy 
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327r3A2m„ Jq^.^ m^zGl + m^G^ ' 
where an upper limit of the integration Gmax is fixed at Gmax('/') = 
2vr^/iamc/^c\/ cos^ Lp + sin^ Lp , where if is polar angle, = nicHzz/mzzlJ^c, 

fizz = _Ua sin^ 9 + fic cos^ 6* is a mass tensor, and an anisotropy parameter 7^ = \fJi^JJia as in 
Ref. |3|. The multi-band torque formula is obtained as 

7!-l sin2g 



where a logarithmic factor a[ff) is given by 



= In 
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em = (sin2^+72cos2e)V2, 0, = e,(^)/7,, e^{9) = {sin' 9+^1 cos' 9^/^ and 6^ = e^{9)/^^ 
3|. Equation ([7]) is reduced to Eq. ([5]) when 7 = 75 = 7a- The multi-band effect is expected 
to take place in oxypnictides 11] and MgB2 [3]. As temperature decreases from down 



to T = 0, a discrepancy between 7^ and 7a becomes appreciable. Theoretically, 7^ increases 
up to ~ 6 while 7^ decreases down to ~ 1 in MgB2 at T = [is], fl^ . 

It has been a long-standing issue that the ambiguity of the parameter 77 is inevitable 
in applying Eqs. and ([7]). The vortex-core contribution to the total energy in high- 
K superconductors is small compared to the magnetic and kinetic energy, and hence is 
neglected in the London model. The energy F of the vortex lattice is replaced with an 



integral from Gmm ~ 27ra~^ with an inter vortex spacing a ~ ^^/(po/B to Gmax ~ 27r^~^ with 
an effective core size ^. The cutoff at Gmax ~ 27r,^~^ to avoid the divergence of Eq. is an 
inherent shortcoming of the London approach. The phenomenological parameter 77 consists 
of a major factor 77' and a correction factor e^"'^ as r] = rj' exp{r]c — 1) |9|. The parameter 77' 
accommodates the uncertainty in defining the core size. The core correction 7/c0o-S/327r^A^ 
is added to the London free energy with an uncertain factor t/^. The London model is not 
beneficial for data analysis because it only gives not a true upper critical field Hc2 but an 
effective upper critical field r)Hc2- It is highly desirable to obtain the Hc2 value by developing 
a new torque theory without containing the 7/ factor. 



ixing the above-mentioned issue on 77 in the preceding London models as 



151 ] introduced a useful cutoff function, and it was followed by Yaouanc 



We succeeded in 
follows. Hao et al. 

et al. {4] to derive the analytical solution of the free energy density for superconductors with 
large Ginzburg-Landau parameter k. The z-component of the flux in the reciprocal lattice 
space is approximated by 
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where = 2{^lGl + ^yGl){l + b^)[l - 26(1 - bf] and Ki{u) is a modified Bessel function 
of the second kind. A reduced magnetic field b{6) = B / Hc2{0) as a function of 9 is given by 



b{9) = {B/H^^') Jsin^ 9 + 7| cos2 9 . 



(10) 



Extending the idea of the multiband London model of the anisotropic superconductors of 
Eq. ([T]) developed by Kogan we express the free energy density F using the local magnetic 
flux hz{Gpq) as 
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where Gr,„ is a discrete reciprocal lattice vector, p and q are lattice indices in a reciprocal 



space 



16j . and hziGpq) is given by 
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using a reduced field b{9). We finally obtain a novel torque formula by angular derivative as 
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where ho{9) = VS%ex/2n'X'^l/\ Vpqi9f = Anb{l + b^)[l - 26(1 - b)^] [cui^xiq - p/2)' + 
p'^/u^x], and uj^xi9) = 2^^ex{9) / \/3'yx£^i9) are functions of 9. Equation f[T5|) is executable by 
numerical differentiation. It is also possible to write down an analytical expression of Eq. (11 3p 
with the aid of the modified Bessel function of the first kind as dvKi{v) / dv = —vKq{v). 

Since our new formula of Eq. (fT3|) does not contain a parameter rj unlike the preceding 
London model Q, it is meaningful to trace the behavior of rj of the London model with 



respect to 7, B, and Hc2- The r] is originally assumed as a factor on the order of unity, and 
subsequently Farrell et al. l^j clarified experimentally as r/ ^ 1.2 ~ 1.5 at temperatures 
near using the magnetization expression M = — 0o/327r^A^ In?7ifc2/-B. The noteworthy 
advantage of using Eq. f[T5]) is that one can directly obtain an upper critical field Hc2 without 
bothering about an indefinite 77 factor. Comparison of two models yields a revisited physical 
interpretation of rj in connection with vortex core and B. 

First, we treat the single band case rather systematically where the condition 7 = 7a = 7^ 
is satisfied. Equation (fT3!) can be reduced to a useful analytical expression as 
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We determine a peak position with respect to angle 6 by using Eq. (fT^ when 7 and 6(6') of 
Eq. f|TO|) are given. We attempt to find the conditions so as to give the same peak angle 6p by 
tuning a parameter rj in the Kogan model while the torque curve is normalized at the peak. 
In the inset of Fig. [H the agreement of the two theories is almost perfect between Op and 90 
degrees while it is not so good between and 6p. In the case of 7 = 7 and i?/iJ^2 = 0-3, we 
find rj = 0.323 at dp = 79.9 degrees. We also carried out the calculations of the torque curves 
under the various different magnetic fields, i.e., B/H^}^^= 0.15, 0.3, 0.6, and 0.8. As being 
remarkable in Fig.[Tl all data in the rj versus /{B /H^}!^) representation are well collapsed into 
a single curve. We consider that this is indicative of the validity of the method to compare 
the two theories by tuning the peak position Op. As 7/(5/7/(^2) increases, 77 increases at 
lower 7/(5/7/^2 ); it forms a peak at 7/(-B/i/j,2) — 10, and it decreases gradually at higher 

Second, we attempt to understand a remarkable profile of Fig. [H Physically the inclusion 
of the upper critical field Hc2 in the Kogan's torque formula comes from the integral interval 
from 2Txa~^ to 271^"^. The Kogan model is based on the presupposition that the vortex core 



cut off always occurs at ^. According to the theory by Hao et al. 15[ and by Yaouanc et 



al. [4], however, the cut-off position is not constant at ^ but varies as .^1,(6). Therefore, the 
major correction parameter rj' essentially comes from the replacement of Hc2 by ri'Hc2 = 
r;'$o/27r,^^ ^ $o/27r{2 = $o/27r,^2 . (^^j^^y^ -^^e approximate 77'"^ as a function of reduced 
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FIG. 1: The parameter rj can be scaled function of -i/{B/W^^) on a single curve, where 7 
is the anisotropy parameter, B is the applied field, and H^^2 is the upper critical field parallel to 
the c-plane. Note that 7-^22 becomes the upper critical field -^^2"^ parallel to the ab plane. The 
inset shows the torque curve of the London model (see the dashed line) of Eq. ([7]) for 7 = 7 and 
i^/i^Jg =0-3, where the phenomenological parameter rj is chosen so as to give the same peak angle 
9p with our model (see solid line) of Eq. (jlSp . 



field h{d) by 



(^v^-^y (1 + 6^) (1-26(1 -6)2) 



r]' \ K ' ^ ' 

where the expression of .^1,(6) was given in Refs. j^, ll5|. In Fig. [21 we show theoretical 77' as 
a function of 1/6. The vortex core becomes that of an isolated single vortex at higher 1/6 
while it compresses at very small 1/6 due to dense packing of vortices. Remarkable is an 
enhancement of 77' at an intermediate 1/6 ~ 3 (see arrow) presumably due to overlapping of 
vortices and resultant spreading of order parameter apart from cores. The core correction 
?7c is essentially the property of an individual vortex, and is not so dependent on 6 and 7. 
This explains qualitatively a particular behavior of 77 in Fig. [H 

Third, the remarkable nature of r/ of Fig. [1] inspires us to make a universal function 
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IG. 2: The qualitative behavior of rj' as a function of 1/6 expected from a theory by Hao et al. 
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15[. The parameter rj' is a major term in the uncertainty parameter rj = r/'exp(7/c — 1) as given 
by Kogan et al. ^, is a vortex core size, and b is the reduced field. 



77(7, B, Hc2) for substituting into the single-band Kogan model of Eq. (jS]). Consequently, the 
Kogan model is also able to give the upper critical field Hc2 without assuming an indefinite 
factor of rj. Details of such an application will be published elsewhere. 

Finally, we briefly discuss the generalized cases where the condition 7^ 7^ 7^ is met in both 
our novel model as well as the Kogan model. In Fig. [3l we show a torque curve obtained by 
carrying out the lattice sum in the reciprocal lattice over + < 10^. We set parameters 
(7a = 2.2, 7^ = 3), (7a = 2, 7g = 5), (7A = 1.7, 75 = 5.3), and (7A = 1.1, 75 = 6). We reveal 
that summation up to p^ + q^ < 10^ is enough to achieve a good convergence while the torque 
curves constructed by using the fundamental spots {p^ + < 1^) is not sufficient to have a 
good convergence. It is numerically confirmed that the condition + q^ < 10^ employed in 
Fig. Ogives actually the same results as those obtained by summing up to p^ + q^ < 50^. 

The phenomeno logical parameter 77 appeared in the Kogan model (see Eqs. ([5]) and (I7j)) 
represents the uncertainty in defining the size of vortex core. In Fig. HI we show the several 
torque curves under the different combinations of two anisotropy parameters 7^, 7a, -?/j.2/-B, 
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FIG. 3: The multi-band curves are obtained by Eq. ()13p for several different combinations of 7a 
and 7^ (see text). The magnetic field is fixed as HI2/B = 3.38 (see Ref. The torque r of 

Eq. (fT3|) is given in units of ^oBV/GAir'^X'^, B is the applied field, H^^ is the upper critical field 
parallel to the c-axis, and p and q are lattice indices in a reciprocal space. 



and 77. The parameters employed here are the same with those in Ref. [3|, but the condition 
of rjH^^^/B = 3.38 is replaced by H^}^!^ — 3-38. We regard 77 as a fitting parameter so as 
to give the same angle at the torque maximum for the case of (7A = 2.2, 7^ = 3) and at 
the torque minimum for other cases both for Figs. [3]and|l] (see arrows). The variation of r] 
depending on the combination of 7^ and 7a may be related with the two competing origins 
in specifying 77 as ?7'exp(?7c — 1) (see above). Further systematic studies of r/ as a function 
of 7a and 7^ are of interest as a generalized case of Fig. [H 

In conclusion, we successfully elaborate the torque formula of anisotropic superconductor 
without containing a phenomenological parameter 77. Unlike the preceding London model 



3[, we are able to estimate the true upper critical field Hc2 by analyzing the torque curve. 
The comparison of our new theory with the London model in the case of 7a = 75 has unveiled 
that the 77 can be scaled very nicely as a function of / {B / H^}!^) . A possible interpretation is 
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FIG. 4: The torque curves using the London model Eq. ([7]) for several different combinations of 
7a and 7^ in units of ^qBV/QAtt'^X^ (see Figs. 1, 2, 3 of Ref. ^) for comparisons). The magnetic 
field is fixed as ie^Hl^/B = 100 (e ~ 2.71 ), B is an applied field, and h''^2 is an upper critical 
field parallel to the c-axis. 



given in view of the field dependence of the vortex core size The behavior of rj was also 
investigated in the multi-band picture of 75 7^ 7a, we found that rj changes its value but it 
still remains on the order of unity. 

We thank M. Ichioka for suggesting the idea to explain Fig. 1 by the field dependence 
of the vortex core size. We are also benefitted by discussions with M. Kato, S. Noguchi, 
S. Kawamata, M. Machida, T. Koyama, T. Tamegai, T. Nojima, and E. H. Brandt. This 
work was partly supported by a Grant-in- Aid for Scientific Research from the Ministry of 
Education, Culture, Sports, Science and Technology of Japan (Grant No. 19206104). 
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